PROCEEDINGS 

OF SCIENCE 



Wind accretion in symbiotic X-ray binaries 



K. Postnov* , N. Shakura , A. Gonzalez-Galan 6 , E. Kuulkers c , P. Kretschmar f 
S. Larsson d , M.H. Finger e / , A. Kochetkova , G. Lu 8 and L. YungelsorV' 

"Sternberg Astronomical Institute, 13, Universitetskij pr., 1 19992 Moscow, Russia 
h Departamento de Fisica, Ingenieria de Sistemas y Teoria de la Sehal,Universidad de Alicante, 
Apdo. 99, 03080 Alicante, Spain 

c European Space Astronomy Centre (ESAC), P.O. Box 78, 28691, Villanueva de la Canada, Spain 
Department of Astronomy, Stockholm University, SE-106 91 Stockholm, Sweden 

e National Space Science and Technology Center, 320 Sparkman Drive, Huntsville, AL 35805, 
USA 

' Universities Space Research Association, 6767 Old Madison Pike, Suite 450, Huntsville, AL 
35806, USA 

8 School of Physics, Xinjiang University, Urumqi, 830046 China 

h Institute of Astronomy RAS, Moscow, 48 Pyatnitskaya Str., Moscow, 119017, Russia 



E-mail: kpostnov@gmail . com| 


nikolai . shakura@gmail 


. corr, 


anagonzalez@ua . es 


Erik . Kuulkers@sciops . esa . int 




Peter . Kretschmar@esa . int 




stefan@astro.su.se 




mark . f inger@nasa . gov 


sunny@sai.msu.ru, guolianglv@sina.com lry@inasan.ru 



The properties of wind accretion in symbiotic X-ray binaries (SyXBs) consisting of red-giant 
and magnetized neutron star (NS) are discussed. The spin-up/spin-down torques applied to NS 
are derived based on a hydrodynamic theory of quasi-spherical accretion onto magnetized NSs. 
In this model, a settling subsonic accretion proceeds through a hot shell formed around the NS 
magnetosphere. The accretion rate onto the NS is determined by the ability of the plasma to enter 
the magnetosphere. Due to large Reynolds numbers in the shell, the interaction of the rotating 
magnetosphere with plasma initiates a subsonic turbulence. The convective motions are capable 
of carrying the angular momentum through the shell. We carry out a population synthesis of 
SyXBs in the Galaxy with account for the spin evolution of magnetized NS. The Galactic number 
of SyXBs with bright (M v < 1) low-mass red-giant companion is found to be from ~ 40 to 120, 
and their birthrate is ~ 5 x 10~ 5 — 10~ 4 yr _1 . According to our model, among known SyXBs, Set 
X-l and IRXS J180431. 1-273932 are wind-fed accretors. GX 1+4 lies in the transition from the 
wind-fed SyXBs to SyXBs in which the giants overflow their Roche lobe. The model successfully 
reproduces very long NS spins (such as in IGR J16358-4724 and 4U 1954+31) without the need 
to invoke very strong magnetic fields. 
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Table 1: Parameters of observed symbiotic X-ray binaries 



1. Introduction 



Symbiotic X-ray binaries (SyXB) form a growing subclass of X-ray binary stars in which 
accretion onto a magnetized neutron star (NS) occurs from stellar wind of the companion red giant 
that does not fill its Roche lobe [jl]]. All these systems are hard X-ray emitters [Q] and are low-mass 
X-ray binaries, with the possible exception of Set X-l [Q], for which a high-mass solution also is 
viable. One of these systems, IGR Jl 6358-4726 may harbour a magnetar Currently about ten 
sources are classified as SyXBs or candidate systems (see Table [I]). 

Remarkable features of SyXBs compared to other low-mass X-ray binaries are their long or- 
bital periods and long NS spin periods. In the catalog of 187 low-mass X-ray binaries [§], the 
orbital periods of 4U 1700+24 and GX 1+4 are the longest in all orbital periods observed. With 
18300 s spin period, 4U 1954+31 is the slowest accretion powered NS [|l|]. Therefore, a detailed 
investigation of SyXB is of immense importance for comprehension of evolution of NS in LMXB. 



2. Quasi-spherical wind accretion 

The cold stellar wind of the red-giant companion is gravitationally captured by NS at the 
characteristic distance Rq = 2GM / (V 2 + v 2 orb ) 2 (the Bondi radius). Here V w is the stellar wind 
velocity at the NS orbital distance, v or j, is the NS orbital velocity. We consider the case of ineffective 
cooling of plasma behind the shock, so the bow shock is formed at the characteristic radius ~ Rq. 
Behind the shock front, the matter falls with a subsonic velocity toward the NS magnetosphere 
located at a distance Ra (the Alfven radius) which is determined by the pressure balance in the 
accreting matter and the magnetic field. To within a numerical factor of the order of a few, which 
depends on the magnetosphere structure, Ra ~ [/I 2 / (M\/2GM)] 2 ^ 7 . 

The matter carries the specific angular momentum j w = k w Q.BR 2 G (here — 1 < k w < 0.4 is the 



numerical coefficient [|17|], D.b is the binary orbital angular frequency). If j w > ]k{Ra) = VGMRa, 
an accretion disc is formed around the NS magnetosphere. In the opposite case, the accretion pro- 
ceeds quasi-spherically. The interaction of accreting matter with the NS magnetosphere ultimately 



leads to the NS spin-up or spin-down - the phenomenon that has been intensively studied [18]. 

The disc accretion on NS magnetospheres is well studied (see JT9[ ] for the discussion) and can 
explain many properties of accreting X-ray pulsars. However, the analysis of spin-down torques in 
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a prototypical SyXB GX 1+4 [20, 21 ] shows that the disc accretion can hardly explain the observed 
correlations. 

The quasi-spherical accretion onto slowly rotating magnetized NSs has recently been revised 



in [|22p . In contrast to Bondi accretion, if the cooling is ineffective, the matter can settle onto the NS 
magnetosphere with a subsonic velocity through a quasi-spherical hot shell. The structure of this 
shell can be found from the hydrostatic equilibrium: p(R) ~ /?~ 3 / 2 , P(R) ~ R~ 5 ^ 2 , T(R) ~ R~ l . 
The hot plasma near the shell bottom can be marginally stable against interchange instabilitites 
mediating plasma entering the magnetosphere, and for accretion to become possible the plasma 
cooling (e.g. via Compton scattering of X-ray photons) may be required [p3|]. However, the very 
fact of observations of slow X-ray pulsars accreting from stellar wind and experiencing spin-down 
episodes suggests that plasma finds ways to enter the magnetosphere. 

In our model the accretion rate onto the NS, which determines its X-ray luminosity, is al- 
ways smaller than the Bondi-Hoyle accretion rate potentially available by the neutron star Mq ~ 
p w R G jv\; the extra matter gravitationally captured by the NS should bend around the hot shell and 
carry away the angular momentum in the accretion wake. 



As shown in [22], the Reynolds number in the hot shell is very high so the interaction of the 
rotating magnetosphere with the shell initiates subsonic turbulence. At the bottom of the shell 
the energy is released due to interaction of the shell with the magnetosphere and turbulence, so 
powerful convective motions begins and establish nearly iso-angular-momentum rotation law in 
the shell: co(R) ~ R 2 . The NS will spin-up or spin-down depending on the specific angular 
momentum of matter at the magnetospheric boundary and the angular velocity of the NS rotation 
00* = 2%jP s . The spin-up/spin-down equation of the quasi-spherically accreting neutron star can 
be written in the form: 

loo* = ZM(o B R 2 G - (Z - z)MR 2 A (0* , (2. 1) 

where Z and < z < 1 are some dimensionless coefficients. Clearly, for spin-down episodes to 
become possible the condition Z > z must be satisfied. Analysis of torque reversals observed in 
wind-accreting X-ray pulsars GX 1+4, Vela X-l, GX 302-1 and Cen X-3 shows that in the case of 



quasi-spherical accretion Z > 1 in these pulsars []22|]. The above equation for the NS spin evolution 



allows us to estimate the NS magnetic moment of X-ray pulsars in which oo* can be measured as 
a function of the X-ray flux near the torque reversal point (i.e. where oo* = 0). Remarkably, the 
estimate of the NS magnetic moment /I is independent of the actual value of the X-ray luminosity, 
only M(doo* / dM) near the torque reversal point matters. The magnetic field of NS in Vela X-l 
and GX 301-2 estimated in this way is consistent with the value inferred from the cyclotron line 
measurements. 

The very slow pulsar spins in SyXBs {oo*Ra < 00k(Ra)) make it hard to establish the pro- 
peller regime to eject matter with parabolic velocities from the magnetosphere during spin-down 
episodes. The important ratio of viscose tensions to the gas pressure at the magnetospheric bound- 
ary is proportional to (oo* — oo(Ra))/oOk(Ra), with the maximum value for the non-rotating enve- 
lope with oo(Ra) = 0. For very slowly rotating pulsars this ratio is very small, i.e. only large-scale 
convective motions with the characteristic hierarchy of eddies scaled with radius can be estab- 
lished in the shell. When oo* > 00k{Ra), the super-sonic propeller regime must set in. In that case 
the maximum possible spin-down torque is ~ /J. 2 /R A - It i s not excluded that a hot envelope will 
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Figure 1: The spin evolution of an accreting magne- 
tized NS in a binary system 
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Figure 2: X-ray luminosities and pulsar spin periods 
in SyXBs with bright red-giants for different models 
(see [^9[ for more detail). Open circles show sources 
from Table 1 . 



persist here, too, and bring away the angular momentum from the rotating magnetosphere. If the 
characteristic gas cooling time in the envelope is short, one can expect the formation of a 'storag- 
ing' thin Keplerian disc [25]. There is no accretion of matter through such a disc, it only serves to 
remove the angular momentum from the magnetosphere. 

At high Mq, thermal instability can appear in hot plasma behind the bow shock front. The 
medium can be stratified into cold dense and hot rarefied regions. Cold blobs move in the hot 
{06T ~ GM/R) plasma. As they move deeper into the hot shell, their velocity can exceed the local 
sound velocity, and blobs get gradually destroyed. Crossing the line of sight, these cold blobs would 
produce an X-ray flux variability with increased hardness ratio. At high X-ray luminosities (high 
M) processes of Compton heating and cooling near and above the Alfven radius become important. 
Below some characteristic distance the Compton cooling by X-rays dominates, and above it the 
Compton heating prevails. This radius can be determined from the condition T(R X ) ~ Tx, where 
Tx is the effective temperature of X-ray photons. The latter depends on the spectral shape (for 
example, T x = 1 /AT for spectrum with an exponential cutoff ~ e~ E l kT , T x = T for the Planckian 
spectrum). If the Compton cooling time is shorter than the fall time of matter, non-stationary 
accretion regime sets in p6|]. Real wind-accreting systems demonstrate complex quasi-stationary 
behavior with dips, outbursts, etc. These features can be interpreted using the proposed model of 
quasi-spherical accretion. 



3. Population synthesis of SyXB 



For the simulation of binary evolution, we use the rapid binary star evolution code BSE [27] 
with updates by [28]. The code is appended by an algorithm for the careful treatment of spin 
evolution of a magnetized NS and the NS magnetic field decay, as described in [29]. The impor- 
tant parameter for the SyXB evolution is the stellar wind velocity from the red giant, which we 
parametrize as v w = ^ p(2GM /R). The parameter j8 is set to 1. 

An example of the spin evolution of a magnetized NS in SyXB is presented in Fig. [jj The 
parameter for the specific angular momentum of captured stellar wind is k w = 1/4. We start with 
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a binary consisting of a 1.4 M Q NS with the surface magnetic field 5 x 10 12 G, the spin period 
P s = 0.01 s, and a 1.3 M main-sequence star. The initial orbital period of the system is Pg = 400 d. 
Before point 'A', the secondary is still a MS-star and Mq onto the NS is below 1O~ 15 M yr -1 , so 
the NS evolves like an isolated object. The evolution of P s is determined by the spin-down torque 
calculated by the magneto-dipole formula. At point 'A the secondary enters the Hertzsprung gap 
and starts to lose mass via stellar wind. The Bondi radius is Rq = 1-3 x 10 13 cm, relations > 
Ra> R c = (GM/ft)* 2 ) 1 / 3 and j w < Jk(Ra) hold, so the NS enters the super-sonic propeller stage of 
quasi-spherical accretion. The NS spin period (and hence the corotation radius R c ) increases, and 
at point 'B' Ra = R c - The NS turns out at the quasi-spherical accretor state from stellar wind, and 
the evolution of P s is determined by Eq. 2.1. The SyXB stage starts. At point 'C, the secondary 
evolves to the first giant branch (FGB) stage. The NS is still in the accretor state. With the ascend 
along FGB, the stellar wind velocity of the giant v w decreases so that j w ~ v w 4 increases. At point 
'D', j w > Jk(Ra), and a wind-fed accretion disc forms around the NS. At point 'E' the secondary 
overfills its Roche lobe, the mass accretion rate immediately increases and exceeds the critical 
Eddington rate. The NS becomes a super-accretor. At point 'F', the secondary leaves FGB and 
evolves to the core helium burning stage. The mass-loss rate rapidly decreases, but the velocity 
of the stellar wind is low so that j w > ]k{Ra) and the NS becomes a propeller with a wind-fed 
accretion disc. At point 'G', the secondary evolves to early asymptotic giant branch (E-AGB). At 
point 'H' Ra < R c and the NS becomes a wind-fed accretor again. For a short time the system 
again appears as a SyXB. Due to the high mass accretion rate, the spin period P s rapidly decreases. 
At point T the secondary overfills its Roche lobe again, so that the supercritical disc is formed. At 
point 'J' the secondary evolves into thermally pulsing AGB-star. At point 'K' the secondary leaves 
the AGB to become a hot white dwarf. To summarize, between points 'B' and 'E' and from 'H' 
to 'K' the considered system appears as a SyXB. Note that the torque reversal from spin-down to 
spin-up occurred close to point 'C when the mass loss rate from the red giant gradually increased, 



in accordance with Eq. 2.1 



4. Results and conclusion 



We have used the treatment of quasi-spherical wind accretion considered in [ 22 ] to calculate 
the properties of a population of SyXB in the Galaxy. The resulting population of SyXBs on the 
P s — L x plane is shown in Fig. [2[ in which known SyXBs from Table [j] are also plotted by open 
circles. Several population synthesis models (with different kick velocities of the nascent NS and 
stellar wind parameters) were calculated (see j j29| ] for more detail). The main results are: 

1. The galactic birthrate of all SyXBs is ~ (4 — 9) x 10~ 5 yr _1 . The total number of galactic 
SyXBs is between ~ 5000 and ~ 10000. The birthrate of SyXB mostly depends on the kick velocity 
dispersion imparted to nascent NS in core-collapse supernovae. However, this parameter affects 
the population of SyXB by a factor of about 2. 

2. 70% — 80% of NSs in SyXB are formed in core-collapse supernovae, 10% — 20% are 
formed in electron-capture supernovae. Less than ~ 1% of SyXBs experience super-Eddington 
accretion. 

3. The galactic number of SyXBs with bright (the absolute stellar magnitude M v < 1) low-mass 
red-giant companion are from ~ 40 to 120, and their birthrate is ~ (5 — 8) x 10~ 5 yr _1 . According 
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to our model, among known SyXBs, Set X-l and IRXS J180431. 1-273932 are wind-fed accretors. 
GX 1+4 lies in the transition from the wind-fed SyXBs to SyXBs in which the giants overflow their 
Roche lobe. The model successfully reproduces very long NS spins (such as in IGR Jl 6358-4724 
and 4U 1954+31) without the need to invoke very strong magnetic fields. 
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